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ABSTRACT 

We improve the bound state approach of the Skyrme model apphed to the heavy 
baryons by adopting a static heavy meson picture where the sohton moves around the fixed 
heavy meson. This allows to take into account the center of mass corrections in a more 
consistent way. The bound state masses so obtained are comparable to the experimentally 
observed Ac and A* masses. A loosely bound state of a soliton with an antiflavored heavy 
meson is found, which leaves a possibility of the nonstrange pentaquark baryon(s). 
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1 Introduction 



In the Skyrme model, heavy baryons containing a single charm or bottom quark can be described 
by bound states of soliton and the heavy meson of the corresponding heavy flavor It was 

first developed to describe strange hyperons [|lO| and was extensively applied to heavy baryons by 
Rho, Riska and Scoccola for which the SU(3) symmetric Lagrangian was strongly modified to 
incorporate large flavor-symmetry breaking. Although qualitatively successful, because of the lack 
of the heavy quark symmetry in their Lagrangian, the resulting heavy baryon spectra could not be 
consistent with the required symmetry. Later, Jenkins, Manohar and Wise have resolved this 
problem by adopting a proper heavy- meson effective Lagrangian |11| that incorporates both the chiral 
symmetry and the heavy quark symmetry explicitly. Based on the model Lagrangian of Ref. [ 11 1 , we 
have carefully examined this picture to study the heavy baryon states with higher spin and also 
to study nonstrange pentaquark exotic baryons jl^Pl. 

These works started out in the large Nc (number of color) and large mq (heavy quark mass) limit 
where the soliton and the heavy mesons are infinitely heavy so that they sit on the top of each other. 
Such a simple picture is not good enough to describe quantitatively the real heavy baryons, especially 
the excited states. In a series of papers [15, 1^, we have tried to improve the model by incorporating 
the finite mass effects of the heavy mesons. We adopted an effective Lagrangian constructed 
for the heavy mesons of finite mass and solved exactly the equations of motion of the heavy mesons 
derived from it. The results are comparable to the quark model calculations in the charm and/or 
bottom sector. 

However, we have worked in the conventional soliton- fixed frame by taking infinite Nc limit but 
with finite niQ. It raises the conceptual problem of 'the tail wagging the dog' [Q]. Actually such a 
problem has been present even in the bound state approach to the strange hyperons. Since the mass 
of the kaon is about one half of the soliton mass, a large center of mass correction is expected. In 
case of the bound state approach to heavy baryons, the situation is worse; the heavy meson is much 
heavier than the soliton — D (B) meson masses are around 2 GeV (5 GeV) whereas the soliton mass 
is about 1 GeV. Therefore it is obviously questionable that the former moves around the latter. A 
phenomenological consequence is that the heavy quark flavor symmetry appears broken much more 
than expected |16|. As a crude way of estimating the center of mass correction, the reduced mass of 
the heavy-meson-soliton system was used instead of the heavy meson mass [16, 18|, which results in a 
quite different spectrum. 

In the real world, however, since the heavy mesons are much heavier than the soliton, the soliton- 
fixed frame may not be a good starting point if one wants to treat any center of mass corrections 
in a perturbative way. It would be more reasonable to start from the opposite picture where the 
soliton moves in the heavy-meson-fixed frame, i.e., infinite mq limit but with finite Nc- Work along 
this direction was first developed in Ref. where the interaction between the heavy meson and 
the soliton was approximated to that of the simple harmonic oscillator. The effect of the kinetic 
motion to the ground state energy is included through the zero point energy. In Ref. ||l^, the same 
harmonic approximation is adopted to describe the excited states Ag. However, the crudeness of the 
approximate potential leads to the Ag states lying high above the experimental observation. 

In this paper, we improve the harmonic approximation of Refs. [|,|l9| by obtaining the equation 
of motion for the soliton in a more realistic way. In the conventional soliton-fixed picture, the equa- 
tions of motion for the heavy mesons have been easily derived from the given heavy meson effective 
Lagrangian after substituting the stationary hedgehog ansatz for the soliton configuration centered 
at the origin. On the contrary, in the picture of fixed heavy meson, we already have a known wave- 
function of the heavy meson which should be localized at the origin. So the static hedgehog ansatz is 
necessarily modified to generate the kinetic motion of the soliton. For this purpose, we introduce the 
collective coordinates representing the translational motion of the soliton. The classical Lagrangian of 

^For an extension of this idea to the NJL sohton model, see, for example, Ref. Jl^ . 



the collective coordinates is then obtained when the given effective Lagrangian density is integrated 
over the space with the heavy meson fields localized at the origin and the chiral field of the soliton 
configuration. Canonical quantization, then, leads to a nonrelativistic Schrodinger equation for the 
motion of the soliton interacting with the heavy meson sitting at the origin. This process is described 
in Sec. 2. 

In Sec. 3, we compare the resulting equation of motion with that of the harmonic oscillator 
approximation. The Schrodinger equation is solved numerically to obtain the bound states. The 
model yields a few bound states of soliton and heavy mesons and the results are compared with the 
available experimental masses of Ac and A* states. We also investigate the bound state(s) of the 
soliton to the antiflavored heavy mesons, which corresponds to the pentaquark {Qq"^) exotic baryons. 
Section 4 contains a summary and conclusion and we provide detailed expressions in Appendices. 

2 Soliton motion in the heavy meson fixed frame 

We start from the heavy meson effective Lagrangian of leading order in l/mg expansion ||ll|, which 
reads 

£ = Cm + iVf,TT[H{d>' + F^^)^] + ggTiiH^'^^^A^H), (2.1a) 

where Cm is an SU(2) x SU(2) chiral Lagrangian governing the dynamics of the Goldstone bosons. 
For a simplicity, we will take the Skyrme model Lagrangian for Cm', 

Cm = ^Tvid^U^d^'U) + ^Trp^d^U, U^d,U]\ (2.1b) 

with the pion field represented by ?7 = exp(ir • t^/ J-k)- The vector and axial- vector fields and 
are defined as 

= ^(e^s^e + e^^e^), (2.1c) 

\ = l{i^d^i-id^£}), (2.1d) 

in terms of ,^ = exp(iT • 7r/2/^). The heavy meson fields are combined into a 4 x 4 matrix field H{x) 
as 

H=^-^{P;i^-Pl^) and H = joHho, (2.1e) 

with the help of the standard Dirac matrices 7^ and 75. Here, P and P* are the heavy pseudoscalar 
and vector mesons moving with 4- velocity v^, respectively. This Lagrangian contains three physical 
constants, the pion decay constant fj^, the Skyrme parameter e and the coupling constant gq of heavy 
mesons. We will take them as free parameters varying around their empirical values (/tt ~ 93 MeV 
and e ~ 4.8) together with a nonrelativistic quark model prediction {gg = — |). One can easily check 
the invariance of the Lagrangian ( |2.1[) under the SU(2)^ x SU(2)j:j chiral transformation 

U^LUR^ {i^Lih^ = hiB)), 

H ^ Hh). ^^-"^ 

In the conventional bound state approaches, the static soliton configuration centered at the origin 
?7o(t") = exp[zr • rF{r)\ is substituted into the Lagrangian and the resulting equations of motion for 
the heavy-fiavored mesons (or kaons) moving in the static potentials arising from the soliton are solved 
to find the bound states. This 'soliton-fixed' picture supported from the large Nc arguments could 
be a good starting point for the strange baryons, which could be improved by proper center of mass 
corrections later. In case of the bound state approach to heavy baryons, however, the situation is quite 



different. The observed heavy mesons D and D* or B and B* are definitely much heavier than the 
soHton. In this case, more reasonable starting point would be the heavy-meson-fixed picture, where 
the soliton moves around the heavy meson fixed at the origin and the corrections from the finite heavy 
meson mass could be included perturbatively. 

In the heavy-meson-fixed frame, the spatial part of the heavy meson wavefunction is already 
determined by a well- localized function at origin such as 6^{r). To generate the motion of the soliton, 
we substitute an ansatz ||2^ 

U{r,t) = Uo{r-x{t)), (2.2) 

into the Lagrangian for the soliton configuration with its center at x{t). Carrying out the integration 
over r, we obtain the classical Lagrangian for the collective variables x{t) asQ 

L= f d?rC = Ms + \MsU^ - M^lit . (/^ X a;) - Vi, (2.3a) 

J X 

where 

Vi = 2gQ[ai{x)Sih ■ Ih + a2{x)Si>h -xlh- x], (2.3b) 

and u = dx/dt is the velocity of the soliton and Mg the soliton mass. (From now on, we will drop 
the trivial soliton mass term from the Lagrangian.) S^h and 1^ are the light quark spin and isospin 
operator of the heavy mesons, respectively. The radial functionals v{x), ai{x) and a2{x) are defined 
as 

sin^(F/2) SinF dF sinF 

v[x) = , ai[x) = and a2[x) = — . (2.3cj 

X X dx X 

The classical Hamiltonian is obtained by taking the Legendre transformation of the Lagrangian 



H,, = ^^^^p' + ^Ih-L + 2v^^ + Vi{x), (2.4) 

with the conjugate momentum p defined as 

dL 2v 

p = — = MsU ilh>ix). (2.5) 

OU X 

Quantization of this collective motion leads us to the Schrodinger equation in the form of Q 

The wavefunction ^'(a;) depends not only on the collective coordinates x but also on the light-quark 
spin sgh and the isospin ih of the heavy mesons. 

Since the Hamiltonian commutes only with the 'light-quark grand spin' operator defined by Ki = 
L + Seh + ^h, the eigenstates can be classified by the corresponding quantum numbers (fc^, m) and 
the wavefunction can be written as 

^{x) = ^ iji{x) {x\ki, m)i, (2.7) 

i 

where the summation i runs over the possible ways of combining the grand spin eigenstates {x\ki,m) 
from orbital angular momentum basis Y£^mi{x), light-quark spin basis \s£h = and the isospin 

basis \ih = i^) of the heavy mesons. 

For a given k£{^ 0), there are four possible eigenstates (see Appendix A) and, because of their 
diff'erent parity, only two states are mixed for a given parity. The wavefunction, then, can be written 
explicitly as 

"^{x) = ipi{x) {x\ki,m)i + 'il^2{x) {x\ke,m)2, (2.8) 
*We treat the soliton motion nonrelativistically. 



for the states of parity tt = —(—1)'^^, and 

"^{x) = tpsi^x) {x\ke, m)3 + ^4(2;) {x\ke, 171)4, 



(2.9) 



for the states of parity vr = +(—1)'^^ including the intrinsic parity of heavy mesons. 

Then, the Schrodinger equation ( ^.61) is reduced to coupled ordinary differential equations for the 
radial functions ipi{x) as 



1 



2M, 



cP 2_d_ 

dx'^ X dx 



or 



with the 2x2 unit matrix I2, where Vqs is in the form of 2 x 2 matrix 

Vesix) -- 

and the corresponding radial function ^/>(x) is 

tjj{x) 

The explicit expressions for the potentials Vij can be found in Appendix A. 
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lp2{x) 



or 



V33 V34 
V43 V44 



^3(2;) 



(2.10a) 



(2.10b) 



(2.10c) 



3 Bound states of soliton 



Before solving the Schrodinger equation numerically, let us compare Eq. (pTO| ) with the har monic 



approximation of Refs. |^, IS]. By making use of the fact that the soliton function F(x) behaves near 
its center as 

F{x) = vr + F'{0)x + lF"'{0)x^ + O(x^), (3.1) 
we can expand the radial functionals of Eq. ( ^.3c| ) as 

v{x) = -- i[F'(0)]2x + O(x3), 

X 

a^{x) = -F'{0) - i{F'"(0) - [F'(0)]3}x2 + 0{x^), (3-2) 

a2ix) = 2F'(0) + l{iF"'{0) - [F'{0)]^}x^ + ©(x^). 

In the extreme limit of Nc ^ 00, the soliton is also infinitely heavy so that it sits on the top 
of the heavy mesons. The binding energy of the system is the value of the potentials at the origin. 
Neglecting the centrifugal energy coming from the kinetic motion of the soliton, we obtain 

V (0) = _Jq£^ ( 2Vf3_ -A^hike + 1) \ 
'^^^ 2{2ke + -Wkeike + 1) 2ke - 1 J' 

for TT = —(—1)'^^ states and vr = +(—1)'^'^ states, respectively. (In order to distinguish the potentials 
obtained with no centrifugal energy, we use the tilde on V.) Diagonalization of those potential ma- 
trices yields the results of Refs. |1^, 13]; that is, degenerate heavy-meson-soliton bound states with 



binding energy — |gQF'(0) and antiflavored-heavy-meson-soliton bound states with binding energy 
— ^5q-F'(0). These degeneracies are of course the artifacts of ignoring the kinetic motion, especially, 
the centrifugal energy. 



Table 1: Vq and k evaluated with conventional Skyrme model parameters. Energies are given in MeV 
unit and k in MeV^ unit. 



Set 




e 


9Q 


Ms 


Vo 


K 


El, 




(I) 


64.5 


5.45 


-0.75 


867 


-791 


(575)^ 


-88 


468 


(11) 


64.5 


5.45 


-1.00 


867 


-1055 


(632)3 


-246 


540 


(III) 


93.0 


4.82 


-0.75 


1105 


-1009 


(733)3 


-113 


597 


(IV) 


93.0 


4.82 


-1.00 


1105 


-1345 


(807)3 


-311 


690 


Expt. 


93.0 






867 


-1095 


(437)3 


-630 


310 



In case of kj = 0"*", we have only one state {x\ki = 0,0)3. Although this state is made of ^ = 1 
angular momentum basis, the associated singularity in the centrifugal energy is canceled by the same 
one of v{x) in the potential, and the resulting effective potential is finite at the origin. Consequently, 
it provides the lowest energy bound state, if any. If we take the harmonic approximation of Eq. (p.2|), 
the potential energy can be written as 

^33(^) -Vq + \kx\ (3.4a) 

where 

Fo = -i5Qi^'(0), (3.4b) 
^ = 9Q{\\F'm'-\F"'m. (3.4c) 

In Table 2, we list the values of Vq and k obtained by using Eq. ( |3.4a| ) with conventional parameters 
of the Skyrme model. 

In this harmonic approximation, the eigenstates can be easily obtained as 



^n = l^o + y^(n + i), (3.5) 
where n is an even integer number. The lowest energy state of n=0 is the ground state which has a 



negative energy, which corresponds to the binding energy of the Aq baryon. From Eq. (3.5), we can 
write the mass formula for Aq(^^) baryons as 



= niN + mu + Vq + ^ ^ [n + fj , (3.6) 

where we have included the mass of the heavy mesons mn and the nucleon mass (soliton mass 
plus the energy from collective rotation) . The kinetic motion of the soliton with finite mass is included 
through the zero point energy of |-\/k;/Ms. The artificial degeneracy in the model of Nc,mQ — > 00 is 
thus removed, while only the ground state is bound. 



In Ref. 1 19], this harmonic approximation of the potential was straightforwardly extended to A;^ > 
states by simply adding the centrifugal energy. The tensor couplings of V12 (V34) between the \ki,m)i 
and \k£,m)2 states (|A;^,m)3 and \ki,m)4 states) was not taken into account in this approximation. 
Then, the potential gives rise to an infinite tower of eigenenergies of Eq. ( |3.5| ) for any integer n. The 
first excited kJ = 1~ states which is ^/kJWs above the ground state could be interpreted as Aq(| ) 
and A^(|") states. In order to fit the observed masses of Ac (2285 MeV), A*(i~) (2593 MeV) and 
A*(| ) (2625 MeV) [^,22], we need to have the binding energy E\, and the excitation energy Ef, as 



= Vb + l^f^W's = -630 MeV, 

(3.7) 



Ee = Vk/Ms = 310 MeV, 




Figure 1: Potentials and wavefunctions in the harmonic oscillator approximation with the parameter 
set (I). The solid line is the exact form of V33 and the long dashed line is its harmonic approximation. 
The potential of kj = 1~ state including the centrifugal energy is given by the dot-dashed line. The 
corresponding wavefunctions are also presented. 

which are regarded as the 'experimental' values in Table |l|. As can be seen in Table |l], under the 
harmonic approximation, it is not easy to fit the experimental values by adjusting the parameters of 
the Skyrme model around their empirical values in a reasonable range, which was also pointed out in 



Ref. HI. 



Shown in Fig. 1 is the harmonic approximation (long dashed line) of the potential V33(x) in 
comparison with the exact one (solid line). The dot-dashed curve is the effective potential in the 
harmonic approximation with the centrifugal energy iesi^eS + l)/(2MsX^) with (£eff = 1) included. 
The lowest energy levels of the corresponding harmonic potentials are explicitly presented with their 
wavefunctions in the parameter set (I). First of all, compared with the exact one, the approximate 
potential is too stiff. Because of this crudeness of the harmonic oscillator potential, the zero point 
fluctuation and the excitation energy become too large. This is crucial for the excited states whose 
wavefunctions are spread over a wide range. Consequently, the binding energy of the Ag appears too 
small and the first excited Aq state locates too high above the experimental observations for all the 
parameter sets adopted in Table |^. 

Now, let us investigate the exact solutions. In Fig. 2, we present two kJ = 0+ bound states 
and their wavefunctions obtained with the parameter set (I). Comparing with those of the harmonic 
approximation (Fig. 1), the exact wavefunctions are more expanded so that the eigenstates get more 
binding. In Fig. 3, solid lines are the potential energies Vij{i,j = 3,4) for the kJ = 1~ state. At a first 
glance, there does not appear any attractive part in each potential, with the exception of V44 which 
is however not strong enough to support a bound state. To understand the binding mechanism, we 
diagonalize the potential matrix at each position x: 

S{x)V^sS\x) = diag(y+, V~), (3.8) 

where 



V^{x) = i[(y33 + ^44) ± ^j{Vzz - Vu)^ + Wi^]. (3.9) 
These are presented by dashed lines in Fig. 3. The off-diagonal potential V34 pushes down V44 so 
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Figure 2: Solutions for the kJ = 0"*" hound states with the parameter set (I). Thick sohd line is the 
potential V33 and the solid lines are the wavefunctions of the ground and radially excited states. 

Table 2: Numerical results on the bound states. Energies are given in MeV unit. 



(n, kJ) 


Set I 


Set II 


Set III 


Set IV 


Exp. 


(0, 0+) 


-287 


-461 


-366 


-588 


-610 


(1, 0+) 


-12 


-62 


-15 


-79 




(0, 1-) 


-89 


-196 


-113 


-250 


-320 


(0, 1+)* 


-17 


-54 


-21 


-69 





* bound state of soliton to antiflavored heavy meson 




that it becomes sufficiently attractive to support a bound state. Of course, one should not use one of 
the diagonalized potentials, or V^, to obtain the eigenstates of the equation of motion because 
the similarity transformation matrix S{x) of Eq. ( p.8| ) is local. Instead, one should follow the general 
procedure to solve the coupled equations as described in Appendix B. The resulting bound state is 
shown in the inner box of Fig. 3 with its wave functions ip3{x) and ip4,{x) (solid lines). There, V~ is 
drawn only as a guiding line. 

We summarize our numerical results in Table |2[ We list the bound states of the model with four 
different parameter sets (I)-(IV). In each case, we could find two bound states in the kJ = 0^ channel 
and one in the 1~ channel. This should be contrasted with the results of the soliton- fixed frame of Ref. 
16|, where many bound states appear as we increase the heavy meson mass. The coupling constant 



predicted by the nonrelativistic quark model (|(7q| = |) yields too small binding energies compared 
with the experimental values. By increasing its magnitude to 1, which may represent the role of light 
vector mesons [18|, we could obtain reasonable values. Comparing with the parameter set (I) and (II) 
with U = 64.5 MeV and e = 5.45, set (III) and (IV) of = 93 MeV and e = 4.82 result in more 
strongly bound states. It is mainly due to the different shapes of the potentials and partly due to 
the heavier soliton mass {Mg = 1105 MeV) of the latter which slightly exceeds the required mass 867 
MeV to fit the nucleon mass. In adopting these parameter sets, we account for the Casimir energy 
[23 1 that will reduce down the final soliton mass. 




Figure 3: Solutions for the kj = 1 hound states with parameter set (I). The solid lines are potentials 
Vij (i,j = 3, 4j and the dashed lines are V^. The wavefunctions are given in the inner box with V~ . 

The last line of Table |2| is the result for soliton bound state to the antiflavored heavy meson, which 
corresponds to the pentaquark {Qq'^) exotic baryon in the quark model where q = u,d quarks. It can 
be searched simply by changing the sign of the coupling constant gq |ll3|, We could find only one 
bound state in the kJ = 1"*" channel. Shown in Fig. 4 are the potentials Vij{i,j = 1, 2) for the soliton- 
antiflavored-heavy- meson interactions (solid lines), the diagonalized one as in Fig. 3 (dashed lines) and 
the wavefunctions il^i{x) and V'2(^) of the bound state (inner box). In case of the nonrelativistic quark 
model value of gq, the state is bound too weakly to survive after the l/mq corrections. However, with 
a coupling constant \gq\ ~ 1 which yields reasonable soliton-heavy-meson bound states, the binding 
energy of the bound state of soliton and antiflavored-heavy-meson is somewhat large and it still leaves 
a possibility of loosely bound nonstrange pentaquark baryons. 

4 Summary and Conclusion 

In this paper, we have carried out the bound state approach of the Skyrme model to heavy baryons 
in a more realistic way. We obtain the bound states of the soliton and heavy meson in the fixed 
heavy-meson picture. It improves the harmonic approximation of Refs. by deriving the soliton- 

heavy-meson interacting potential in a more consistent way from a given effective Lagrangian. First 
we introduce the collective coordinates for the motion of the soliton. Then, the classical Lagrangian 
for these coordinates is obtained by carrying out the integral over the space with localized heavy 
meson fields. The nonrelativistic Scrodinger equation is obtained by quantizing the motion. The 
resulting bound states have larger binding energies than in the harmonic oscillator approximation. 
Apart from the ground state with kj = 0~^, we could find the first excited state with negative parity 
in the kJ = 1~ channel. With the value oi gq ~ 1, which is slightly larger than the nonrelativistic 
quark model prediction, those resulting bound states are comparable to the experimentally observed 
Ac{^'^) and A*(^~,|~) states. 

This work also improves the work of Refs. [16| where the bound state approach was carried out 



in the conventional soliton-fixed picture. In this approach the heavy quark flavor symmetry is rather 
strongly broken as expected. Since the l/mq order terms are not included in the Lagrangian, the 
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Figure 4: Solutions for the kj = 1"*" soliton-antiEavored-heavy-meson bound states with the parameter 
set (I). The solid lines are potentials Vij (i,j = 1,2) and the dashed lines are V^. The wavefunctions 
are given in the inner box with V~ . 



obtained bound states do not depend on the flavor of tlie lieavy mesons at tliat level of approximation. 
If one turns on the I/tuq corrections, this degeneracy will be removed as much as the heavy flavor is 
broken by the l/mq order terms in the Lagrangian. 

We have obtained just those bound states which are eigenstates of the light-quark grand spin 
operator. They could be naively compared with the Aq{^~^) and Ag(i ,| ) baryons. However, in 
order to provide those bound states with the correct quantum numbers such as spin and isospin, we 
have to quantize the zero modes of the bound states under the SU(2) rotation in isospin space. This 
procedure will enable us to check the model by comparing the results with the other observed heavy 
baryon states such as Eg and Eg [^,25]. Work in this direction is in progress and will be reported 
elsewhere 
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A Matrix elements in kf basis 



In this Appendix, we give the explicit expressions for the potential Vij of Eq. ( 2.1[1| ). 
The light-quark grand spin operator is defined as 

Ke = L + Ih + Sm, (A.l) 



where Ih and Sih are the isospin and light-quark spin operators of the heavy mesons and L is the 
orbital angular momentum operator of the soliton. It should be distinguished from the grand spin 



Table 3: Four {x\k£m) eigenstates for a given k^. 



parity vr = -(-1)'^^ 


parity vr = +(-1)'^* 


state 


^^ 


\ 


state 




Ai 


{x\ki,m)i 
{x\ke,m)2 


Aj — 1 


ke + ^ 
ke-2 


{x\ke,m)3 
{x\ki,m)4 


Aj + n 

Aj - 2 


ki + ^ 
ke-2 



operator K that is the sum of the orbital angular momentum and all kinds of the spin or isospin of 
the particles in the model. Note that in the definition of the operator Ki only the angular momentum 
operators of the light degrees of freedom are included. 

The eigenstates of the light-quark grand spin operator can be obtained by linear combinations of 
direct product of the eigenstates of each angular momentum operator in its definition: 



{x\k^ 



m's 



'{ii mi^mtl Ximx){Xi mx^ms \ A;^ m)Yf-m^(a;) \^mt) 



(A.2) 



with the help of the Clebsch-Gordan coefficients [27|. Here, we first combine L and Ih (A and mx 
are the corresponding quantum numbers of the eigenstates of A = L + 1^), and then combine the 
light-quark spin Sih- For a given ki{> 0), there come out four eigenstates as listed in Table ^. In case 
of ki = 0, there can be only two states {x\ki, 171)1^3. In evaluating the total parity of the states as in 
Table ^, the negative intrinsic parity of the heavy mesons is taken into account. 

To obtain the potentials Vij{i, j=l,2 or 3,4) of Eq. ( |2.10| ), we need to evaluate the expectation 
values of the operators, Ih'L, Ih-Sih and {Ih-x){Sih-^), with respect to the eigenstates \ki,m)i{i =1,2 
or 3,4). The first one is simple to evaluate the expectation values, since it can be rewritten as 



Ih L 



-I h 



1/2, 

-(^ + l)/2, 



if A 
if A 



Thus, it is trivial to obtain ■ L)j as 



+ 



(A.3) 



2- 





1 


2 


hj 


3 


4 


1 


kki 





3 


-Uke + 2) 





2 







4 





Uki-'^) 



(A.4) 



The expectation values of the operator I^-Sih can have a simpler form if we obtained the eigenstates 
of the operator in another way; that is, by combining Sih and Ih first (S^ = Ih + Sgh with the 
quantum numbers a and m-o-), and then with L. Let \i, a; ki, m) be such a state. Then, the eigenstates 
|A;£, m)i can be related to the eigenstates |£, ct; kg, m) by the Racah coefficients as 



\h,m)i 
\ki,m)2 



h + l 



\i = k£,a = 0; ki, m) + 



\i = kf ,a = 0; ke, m) + 



\l = ki,a = 1; ki,m). 



+ 1 



]] 2ki + l 
\ke, 111)3 = \i = ki + l,a = l;ki,m), 
\h, 111)4 = \i = ke - l,cr = l;ke,m). 



2kp + 1 



\i = ki,a = l;ke,m) 



By using the identity 



Ih ■ Sih = -^i^e - ■^), 



(A.5) 



(A.6) 



we can obtain i{Ih • Sih)j as 
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(A.7) 



Subtracting ^{Ih • -S'^^), the third operator can be expressed in a tensor operator form of 

Ot = {Ih ■ x){Seh • x) - • Sen). (A.8) 

By rewriting it as 

Ot = 2{Ih ■ x){Ih ■ S,h){Ih ■ A) + lilh ■ Sih), (A.9) 

= —\k£,m)i^2 (^=1)2), we obtain i{OT)j as 



which enables us to use the relations 2(i"^ • x)\k£,m)i 
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hi 


3 
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ki 


^/hiki + 1) 


3 


/c£ + 2 


V^K^^ + 1) 


6{2ke + 1) 


6{2ki + 1) 


6(2A;£ + 1) 


2(2A;^ + 1) 
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^ke{ke + 1) 


ke + 1 
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^Jh{ki + 1) 


h-i 


6{2ke + 1) 


6{2ki + 1) 


2(2fc^ + 1) 


6{2ki + 1) 



(A.IO) 



Finally, we obtain Vij{i, j=l,2 or 3,4) as 

Vu =Vr + Vcent + ^VjL 



2ki + 3 



2 4(2fc^ + l) 

2A;^ + 1 6(2A;£ + 1 

t/22 — Vr + Vcent ViL 



ke 



6{2ke + 1 
ke + 1 



2 

ke + 2 



^33 — Vr + Vcent 

2(2fc^ + l) ^' 

V44 = + Vcent + 



4(2fc£ + 1) ^ 6(2fc^ + 1) 

1 kf + 2 

ViL + -:Vis - .J , ,, Vt 



Vt, 



(A.ll) 



6i2ke + 1) 



z-^^^ + ^^^^- 6(2^ + 1)^^- 
where V- = v^Ms, Kent = £{£ + l)/{2MsX^), Vil = 2v/{Msx), Vis = ^gqiai + ias) and Vt = 2gQ02. 

B Numerical Process 



Here, we briefly describe the numerical process of finding the bound states. Equation (2.10) can be 
solved like any other eigenvalue problem. First, we find the asymptotic solutions near the origin and 
at large x. Except the ki = case where only one radial function (tlji{x) or 7/^3 (x)) is defined, we 
can find two independent solution sets. (See Table |^.) Taking those independent sets of asymptotic 
solutions as initial conditions, the differential equation can be solved numerically in the forward or 
backward direction of x. In this way, we can obtain 8 radial functions (2 independent radial functions 
X 2 independent initial conditions x integrated forward or backward direction of x). Let us distinguish 
those radial functions by using scripts >(or <) and A{ot B). For example, "^j^^ denotes the radial 



Table 4: Two independent asymptotic solutions near the origin and at large x. 







near the 


TT = 

origin 


at large x 


vr = 
near the origin 


at large x 


Sol. 


A 


ipi{x) ~ 

lp2{x) ~ 





'il^2{x) ~ 


^'3(2;) ~ x^^ 

Mx) ~ 


V'4(x) ~ 


Sol. 


B 


iplix) ~ 

lp2{x) ~ 





V'i(x) ~ 

4^2{x) ~ eV^^/x 


^p3{x) ~ 

V'4(x) ~ X^'^ 


tp3{x) ~ 

ip4{x) ~ e^^l^^/x 



function obtained by starting from the asymptotic solution set A and solving the equation of motion 
in backward direction. 

Then, a general solution that one can obtain by integrating the equation of motion in forward 
direction is expressed as 

f \ ^( (^Ai^i,A + oiBi^i,B \ 



(B.ll 



where ^ are arbitrary constants. Similarly, a general backward solution can be obtained as 

( ^nx) \ ^ ( pa^Ia+Pb^Ib \ 



(B.2) 
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\ i^f'Ai^m) 


i^2,B 


yXra) 


^2,A(^m) 


i^l'B 


yXm) 
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-Pb 
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At a proper intermediate position Xm, those wavefunctions should be matched in the way that the 
functions and their derivatives are continuous. This matching condition leads to a linear equation. 



(B.3) 



In order that it supports a nontrivial solution, the determinant of the matrix should vanish. Such 
a secular equation determines the eigenenergies. Once an eigenenergy is found, Eq. (|B.3| ) can be 
solved for the ratios between the coefficients, aA,B and Pa,b- Their absolute values can be fixed by 
normalizing the wavefunction. 

In the special case of ki = 0, there appears only one radial function, 'ipi{x) or ^^{x). The asymptotic 
solution of them is given as the solution set A in Table H. In a similar way described above, one 
obtains the backward and forward solutions, V'i(3) and il^ff^^y respectively. Simply by matching their 
logarithmic derivatives as 



(3) 



(B.4) 



one can find the eigenstates. 
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